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Abstract 

When gas atoms or molecules collide with clean and ordered surfaces, under many circumstances 
the energy-resolved scattering spectra exhibit two clearly distinct features due to direct scattering 
and to trapping in the physisorption well with subsequent desorption. James Clerk Maxwell is 
credited with being the first to describe this situation by invoking the simple assumption that when 
an impinging gas beam is scattered from a surface it can be divided into a part that exchanges 
no energy and specularly reflects and another part that equilibrates or accommodates completely 
and then desorbs with an equilibrium distribution. In this paper a scattering theory is developed, 
using an iterative algorithm and classical mechanics for the collision process, that describes both 
direct scattering and trapping-desorption of the incident beam. The initially trapped fraction 
of particles can be followed as they continue to make further interactions with the surface until 
they are all eventually promoted back into the positive energy continuum and leave the surface 
region. Consequently, this theory allows a rigorous test of the Maxwell assumption and determines 
the conditions under which it is valid. The theory also gives quantitative explanations of recent 
experimental measurements which exhibit both a direct scattering contribution and a trapping- 
desorption fraction in the energy-resolved spectra. 
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I. INTRODUCTION 



Trapping and sticking are important processes that occur in gas-surface interactions. 
Sticking is generally associated with strong chemical bonding to the surface and once bound 
the stuck particle does not desorb. Trapping, on the other hand, is associated with the 
physisorption potential well created by the relatively weak Van der Walls potential and in 
many circumstances trapped gas particles will desorb after a period of residing in the well 
near the surface. 

Maxwell, in his studies of gas-surface interactions, is credited with being the first to invoke 
the simple assumption that a gas impinging on a surface is scattered into two fractions, 
one that reflects specularly and exchanges no energy and the other that equilibrates or 
accommodates completely and desorbs with an equilibrium distribution.- This idea was 
taken up early in the twentieth century by Knudsen who introduced the concept of the term 
"coefficient of thermal accommodation" to measure the efficiency of energy exchange at the 
interface between a gas and a surface and developed a theoretical framework in which to 
describe it.- Since this early work it has become standard to assume in gas-surface collisions 
that the fraction of the incident gas beam that is trapped and subsequently desorbed leaves 
the surface in an equilibrium distribution, i.e., its accommodation coefficient is assumed to 
be unity and its distribution function is the Knudsen flux. This assumption of Maxwell 
appears to be very useful because it appears to explain qualitatively experimental results 
measured under a wide range of different conditions, however, such an assumption has never 
been adequately verified theoretically. 

In fact, it is clear that the first part of the Maxwell assumption, i.e., that the direct scat- 
tering contribution is elastic can hold strictly true only in the case of quantum mechanical 
conditions which implies low energies, small temperatures and small mass ratios. For classi- 
cal scattering, which is the regime of large energies, high temperatures and large mass ratios 
the direct scattering contribution leaves the surface as a distribution over a range of energies 
whose average is typically smaller than the energy of the incident beam, but may be larger 
in the case of high surface temperatures.-'^'^'^'^ However, the second part of the Maxwell 
assumption, i.e., that the trapping-desorption fraction leaves the surface in equilibrium, is 
often still today used to explain experiments.-!^ 

The purpose of this paper is twofold. First, we test the Maxwell assumption of equilib- 



2 



rium for the trapping-desorption fraction using realistic calculations for a simple model of the 
gas-surface potential and determine when such an assumption is valid. Second, we demon- 
strate that a model of the interaction potential that retains the basic elements necessary for 
trapping, when combined with a calculation that contains correct statistical mechanics, can 
explain modern high-precision energy-resolved scattering measurements. The gas projectile 
is taken to be an atom and its interaction potential is taken to be an attractive square well 
with a strongly repulsive surface barrier. The calculations are carried out using classical 
mechanics, which is justified for many systems of interest in rarefied surface dynamics. A 
classical treatment means that the results will describe heavy mass atoms at higher energies 
and surfaces at high temperatures where quantum mechanical effects are not dominant. The 
use of an attractive square well to approximate the slowly-varying Van der Waals potential is 
also reasonable when used with a classical calculation since it gives a good description of the 
two primary effects upon entering the well which are an increase in energy and a refraction 
of the atom to steeper angles toward the surface. The authors have earlier presented results 
using a similar model for one-dimensional scattering. This work extends those earlier results 
to the much more realistic and more complicated case of fully three-dimensional scattering.- 

Once a beam of incoming atoms interacts with the surface, a fraction will be directly 
scattered while the remainder will be trapped in the potential well. Of the trapped fraction, 
some will lose sufficient energy to be actually trapped in the well with negative total energy 
while others, even though they have positive total energy, will scatter at angles sufficiently 
close to grazing that they will be deflected back towards the surface by the attractive part 
of the well. This latter positive energy part of the trapped particles is often called the 
chattering fraction. The trapped portion of the incident beam particles will continue to have 
interactions with the surface and with each subsequent interaction some will receive enough 
energy and will be projected sufficiently close to the surface normal that they can escape, 
while the remainder will continue to be trapped. Eventually, in a closed system, all particles 
initially trapped will ultimately desorb from the surface, although for low temperatures and 
deep potential wells this could take a very large time. 

For the model considered here, through an iterative algorithm, all of the trapped particles 
are followed as they continue to have collisions with the surface, and at each iteration 
the negative energy fraction, the chattering fraction, and that fraction which is desorbed 
is recalculated. In this manner the energy distribution of the slowly diminishing trapped 
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particle fraction as well as the energy and angular distribution of the desorbed (or scattered) 
particles can be followed, and a trapping time can be calculated. By following the initial 
direct scattering and the sum of all the subsequently desorbed particles the approach to 
equilibrium of the trapping-desorption fraction can be monitored. 

What is determined is that for shallow potential wells and large surface temperatures, 
conditions for which little trapping is expected, the desorbed fraction leaves the surface very 
quickly and does not at all resemble an equilibrium Knudsen flux. However, for deeper wells 
and lower temperatures when the majority of the incident particles are trapped, it is found 
that the energy distribution of the trapped fraction rather quickly saturates to a stable 
functional form while the total number of trapped particles slowly diminishes. As a function 
of initial conditions favorable to trapping, e.g., low incident beam energy, deep wells and 
low surface temperatures, the approach of the scattered particles towards an equilibrium 
distribution is followed. We find that in many circumstances, the energy distribution of the 
scattered particles readily approaches equilibrium shape even for trapping times that are 
relatively short. The approach of the angular distribution to equilibrium shape, which is a 
Knudsen cosine distribution independent of azimuthal angle, occurs more slowly and only 
for very large trapping times. Thus this work provides a real prediction for the conditions 
under which the Maxwell assumption of equilibrium for the trapping-desorption fraction can 
be applied with reasonable accuracy. 

However, the mere fact that these calculations can indicate the conditions under which 
the trapping-desorption fraction may appear as a nearly equilibrium distribution is not 
sufficient to demonstrate that such conditions are realistic. In order to be convincing, cal- 
culations with the same potential model should be capable of explaining real experimental 
measurements. To demonstrate this ability, we have chosen to compare our calculations 
with recent high-precision energy-resolved data for the scattering of Ar atoms from a well- 
ordered monolayer of the polymer 1-decanethiol adsorbed on Au(OOl).^ This experiment 
showed that for well-defined beams of Ar incident over a large range of energies and angles 
the scattered distributions could be described by a combination of two features, a direct 
scattering fraction and a trapping-desorption fraction that was nearly in equilibrium with 
the surface. Our calculations provide an excellent description of both contributions of the 
scattered spectra and produce an estimate of the average physisorption well depth of the 
interaction potential. 
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In the remainder of this paper the theory is fully developed together with a description of 
the iteration method in the following Section [III In Sec. IIIH a number of calculated results 
describing the approach of the trapping-desorption fraction towards equilibrium are shown 
and discussed. In Sec. IIV[ calculations are compared with the experimental data of Ref. \^]. 
Conclusions are discussed in Sec. El 

II. THEORY 

A convenient way of describing a surface scattering event is through a differential reflection 
coefficient, written as dR{pf,pi)/dEfdQf which gives the fraction of an incident beam of 
momentum pj that is scattered into the small energy interval and small solid angle in the 
direction of the scattered momentum pf. The differential reflection coefficient obeys the 
unitarity condition which assures that the number of particles scattered equals the number 
incident on the surface 



The interaction potential is specified by a vibrating repulsive wall with a square ph- 
ysisorption well in front of depth D and width b, where the actual length of b plays a role 
in calculating the trapping times but is unimportant for all other calculations as long as 
it is larger than the selvage region containing the vibrational corrugations of the repulsive 
potential. It is assumed that the exchange of energy and momentum with the surface occurs 
only in collisions with the repulsive wall, while all trapped particles collide with a stationary 
wall at the front of the well positioned at z = b that simply reflects specularly. This means 
that for the purpose of calculating the resulting distribution after a collision the differential 
reflection coefficient is calculated using momenta that include the well depth in the normal 
component. For example, a particle that would have momentum outside the well has the 
momentum p^ inside where the two differ in that the energy associated with perpendicular 
motion is increased by the well depth 



where m is the atom mass. Similarly, an atom incident from asymptotically far away with 
polar angle 6i with respect to the surface normal will be refracted inside the well into the 
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angle 

l EiCos^e, + D 
^"^^- = V E, + D • 
The relationship between the differential reffection coefficients inside and outside of the well 

for the fraction which escapes is given by a simple Jacobian that is calculated from Eqs. ([2]) 

and ©. 

Within this model, the trajectory of a given atom consists of successive collisions with the 
surface. The incoming beam first enters the well and then proceeds to have a first collision 
with the repulsive wall that scatters into a distribution of energies and angles dictated by 
the differential reflection coefficient. Some of these scattered particles have sufficient energy 
and small enough polar angles to escape out of the well, this is the direct scattering portion. 
The remaining particles are trapped, they are specularly reflected by the front face of the 
potential well, and then they travel back to the repulsive potential where they suffer a 
second collision. This process repeats multiple times until all of the initially trapped atoms 
eventually escape the confines of the potential well. 

Based on a zeroth order differential reflection coefficient dR^{pf,pi)/dEfdQf which, for 
each collision with the repulsive potential, gives the probability of scattering from momentum 
state Pi top f the total differential reflection coefficient after n such collisions can be written 
schematically as 

dR''{Pf,Pi) ^ dR^{pf,Pi) r dR%pf,pb) dR%p,,p,) 

dEfdQf dEfdQf J ^ ^ dEfdQf dEbdQb 

I f^T^^n dR^jPf^Pb) dR\pb,Pi) 

+ J ^^^^^^ dEfdnj dEA + ■■■ 

, -O ^-R°(P/.P^) dR^-\p,,Pi) 

+ J at.,au, ^^^^^^ ^^^^^^ 

where the intermediate integrations in the higher order terms are carried out only over those 
energies and angles that pertain to particles trapped in the bound states. 

Such a procedure lends itself to an iterative formulation in which the scattered distribution 
remaining in the well after the last collision becomes the source for the next collision. The 
angular and energy space within the well is divided into bins sufficiently small so as to obtain 
good numerical precision and it is necessary to keep track separately of the three different 
types of trajectories, i.e., the trapped particles with negative total energy (the trapping 
fraction), the trapped particles with positive energy (the chattering fraction) and those that 
escape at each iteration (the trapping-desorption fraction). 



An explicit mathematical description of how this is accomplished is as follows. After 
the nth {n > 1) iteration, the differential reflection coefficient inside the potential well 
dR''{p'f,p'i)/dE'fdn'f is: 



r (ii?"-Vp'„ p') dR^ 



dE'jdn'j 



dE'jdQ'j dE'fdQ'j ' ^ ' / 



(5) 



where the upper line of the left hand side of Eq. ([5]), labeled with the conditions Ef>D and 
< 6'jr < gives the intensity scattered into the continuum states after n iterations and 
consists of the contribution that was already in the continuum state after n — 1 iterations 
plus the fraction contributed to the continuum by the nth iteration. The critical angle for 
reflection of particles in the positive energy chattering fraction from the front of the well is 
which is dependent on energy and given by an equation similar to Eq. (IHl). The fraction 
that remains trapped in the well is divided into the sum of the positive energy chattering 
fraction and the negative energy trapped fraction, denoted respectively by the subscripts C 
and T, according to 

^^"(P/,Pi) ^ dRUp'f,pd dRUp'f,pd , . 

dE'jdn'j dE'jdn'j dE'jdn'j ^ ' 

The positive energy chattering and negative energy trapped fractions are further divided as 
follows: 

dRWfM) f. . , dRl-\p'f,p[) dR-.cMM) dR-,^^cMM) \ 1 

dE'^d^i'j y ^'^ dE'fdn'j dE'fdn'f dE'fdn'j j ' 



(7) 



and 



dRUp'f,pd _L ,r( >.n>. /RT-\p'f,Pd C^/g7c(T)(PpP:) ^/g,V) (Pp PD ] 1 

dE'jdn'j Y dE'jdn'f dE'jdn'^ dE'^dn) ]w' 

(8) 

In the above equations the factor [1 — N{p'j\9'f)\ where N{p'j\d'f) is the ratio of normal 
velocity to that of the maximum normal velocity of all bound states, given by 

pcos9 

z 

where P^°-^ is the largest normal momentum component of all the trapped particles and 
pcos6 is the momentum component in the z direction for any other trapped particle. The 
term multiplied by this factor takes account of the fact that the slower particles collide 



less often than the faster particles. The TV" is a normalization coefficient chosen such that 
dR'^{p'pp'j)/dE'j:dQ'j is normalized as in Eq. ([1]). The intermediate differential reflection 
coefficients for the chattering and negative energy trapped fractions at each iteration are 
given by 

dE'dn' - Id L L ^'^^ dE'dn' ^^P^'^^^ dE''dn'' ' 

and 



ciE^rffi^ - X ^^^X ^^^io ^'^^^EpfiT^ ^^^^^^ ' ^ ^ 

where the symbol X can stand for any one of the three possibilities: C for the chattering 
fraction, T for the negative energy trapped fraction or Con for the fraction that goes into 
the continuum. For example, dK^rp^^j^^p'j:, p[) / dEjdfl'j is the intermediate differential reflec- 
tion coefficient giving the probability during the nth iteration that a particle will make a 
transition from the negative energy trapped fraction to the chattering fraction. Finally, the 
contribution to the continuum states in the total differential reflection coefficient of Eq. ([5]) 
coming from the nth iteration is 

dR^ConiP'f^P'i) _ 1 dRwiConM^P'i) , 1 ^^?T(Con) (P/> PO 



dE'fdQ'j A/"" dE'^dQ'^ AT" dE'fdQ'^ 



(12) 



At the end of n iterations the fraction of all incident particles that remain trapped in the 
positive energy chattering states is 

P ff ^2.^ dP^{p'f,p'^ 

ft: J J 

while the fraction trapped with negative total energies is 

^- = 70 "^""fJo "^^fJo ^'^z dE',dn', ■ 

Thus the total trapped fraction after n iterations is 

P'' = P^ + P^ . (15) 
The fraction escaping into the continuum state after n iterations is 

,16, 
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and the unitarity condition assures that the total number of particles is conserved 



P5 + P^ + P5„„ = l. (17) 

The major numerical operation in this procedure are the two volume integrals associated 
with evaluating the intermediate differential reflection coefficients of Eqs. (11 01) and ffTTj) 
and the unitarity summations of Eqs. (fT3!) . (IT^ and (fT6!l . which taken together amounts 
to a 6-dimensional integral. The angular integrations are carried out using Gauss-Legendre 
quadratures and the energy integrals use Gauss-Laguerre quadratures. Because classical 
differential reflection coefficients are positive definite and typically tend to consist of a single 
broad peak or a small number of such peaks in both the energy and angular variables, Gauss 
quadratures are ideally suited for these integrals. 

The only remaining element of this procedure that needs to be specified is the zeroth order 
differential reflection coefficient. There are a number of choices that have been used in the 
past to describe classical mechanical collisions of atoms with vibrating surfaces .-'^^I'^i The 
simplest of these, and the most appropriate for the present calculations, is the differential 
reflection coefficient for an atomic projectile colliding with a surface of discrete scattering 
centers of mass M whose initial momenta are distributed in an equilibrium distribution at 
temperature T5. This is given hy^'^^^^: 

dEfdQf Sn^h^p,,^ KkBTsAEoJ ^ \ 4kBTsAEo J ' ^ ^ 

where AEq = (p/ — pj)^/2M is the recoil energy, piz is the z component of the incident 
momentum, ks is the Boltzman constant, \Tfi\ is the form factor of the scattering center 
which depends on the interaction potential. To lowest order, the amplitude rjj is identified 
as the transition matrix element of the elastic interaction potential extended off the energy 
shelU^, however, for this work we use the value appropriate for hard sphere scattering which 
is a constant. The differential reflection coefficient of Eq. (fT8|) can be obtained from a purely 
classical calculation or from a quantum mechanical formulation in which the classical limit 
is extracted. In the case of a completely classical derivation the constant h is unspecified 
except for its dimensions of action, whereas quantum derivations identify h as Planck's 
constant divided by 27r. 

To obtain the trapping time r a variety of methods can be used, but we have found that 
the most convenient is to first calculate the average speed normal to the surface for the 
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trapped particles after each iteration. The time for that iteration is then determined as that 
required to travel the distance 2b from the repulsive wall to the front of the well and back. 
The average normal speed for the positive energy trapped fraction is 



and for the negative energy trapped fraction it is 



The total trapping time is then given by summing the times at each iteration 



T 



( 1 P" 1 P"\ 

^"Hy—^-^^^—^fA + (21) 



where P° is the fraction of initially trapped particles. The actual definition of trapping time 
used in this work is the time required for the fraction of trapped particles remaining in the 
well to be reduced to one percent of the number of incident atoms. 

The method of calculation of the average trapping time is clearly not unique and we have 
evaluated it several ways using the trapped fraction probabilities as in Eqs. f|T9l) and fl20|) 
above. For example, instead of determining the average speed one can use the root mean 
square normal speed, or find the average time directly by obtaining the average of Ihjvz at 
each iteration. In cases in which the trapping time is relatively long, all of these different 
methods yielded values which were quite similar. For all average trapping times reported 
here the width of the well was taken to be 6 = 3 A and Eq. (|2T1) shows that r scales linearly 
with h. 



III. APPROACH TO EQUILIBRIUM 

In this section we present a number of calculations that demonstrate the approach to- 
wards an equilibrium distribution of the trapping-desorption fraction. The parameters are 
primarily chosen to represent monoenergetic and angularly defined beams of either Ar or 
Ne scattering from a tungsten surface. This leads to a set of guidelines for when one may 
expect the Maxwell assumption to be valid, i.e., for when the trapping-desorption fraction 
approximates an equilibrium distribution. 
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Fig. [T] shows an example calculation of the evolution of the energy distribution as a 
function of number of iterations for the case of argon scattering from a tungsten surface. 
The incident angle is 45°, the incident energy is 1 meV and the well depth is chosen to 
be 80 meV and the surface temperature is 303 K. The dotted curve shows the continuum 
energy distribution after the first iteration, which is the second collision with the surface. 
The trapping fraction is = 0.953, indicating that 95.3% of the incident particles remain 
trapped in the potential well. The dashed and dash-dotted curves show the evolution of 
the continuum scattered distribution after increasing numbers of iterations of 5, 50 and 500. 
After 500 iterations there is still approximately one third of the incident particles trapped. 
After 2124 iterations the trapped fraction drops below the arbitrary threshold of 1% of the 
incident particles. 

Also shown for comparison in Fig. [1] is the Knudsen distribution given by 



It is clearly seen that the total scattered distribution closely approximates the Knudsen flux 
when nearly all the initially trapped particles have been desorbed. 

What is actually plotted in Fig. [1] is the integral over all final angles of the differential 
reflection coefficient, which is essentially the average energy distribution scattered over all 
outgoing angles. However, the energy distribution at any fixed final polar and azimuthal 
angle behaves quite similarly and at the maximum iteration number the dependence of the 
distribution on energy is essentially the same at all angles. This is in agreement with the 
Knudsen fiux which has exactly the same energy dependence at all final angles. 

The average trapping time as calculated using the average normal speed from Eq. (fT9ll 
in order to reach the arbitrary threshold of 1% of the particles still remaining trapped is 
r ^ 4.4 X 10^^ s. If the number of iterations is extended to larger than the = 2124 shown 
in Fig. [1] there is essentially no change in the scattered distribution because the number of 
particles remaining in the well is insignificant. 

For incident energies small compared to the well depth D the final, converged energy 
distribution after a large number of iterations is independent of incident energy and incident 
angle. However, the distributions calculated after only a small number of iterations will vary 
somewhat with the choice of these incident parameters. 

The evolution of the angular distribution for the same Ar/W system as shown in Fig. [1] is 




(22) 
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given in Fig. [2l It is seen that as the number of iterations increases the angular distribution 
as a function of the polar angle 6f gradually approaches the cos6f form of the Knudsen 
flux which is shown as open circles. The calculations of Figs. [T] and [2] very quickly become 
independent of the azimuthal angle even for very small numbers of iterations. This is 
consistent with the behavior of the Knudsen flux which is independent of azimuthal angle. 

For the set of initial conditions chosen in Fig. [2] it is seen that the angular distribution 
does not achieve a Knudsen distribution even at the largest iteration number calculated. 
The shape closely resembles a cosine function but, for example, its value at 0j = is only 
about 88% of that for the corresponding Knudsen cosine. This behavior is typical of many 
of the systems examined, if the well depth is sufficiently large, and the incident energy small 
compared to D and ksTs the system will usually achieve an equilibrium energy distribution 
after large numbers of iterations and it becomes independent of azimuthal angle very quickly, 
but the polar angular distribution is usually very slow to converge to the equilibrium cosine 
distribution as seen in Fig. [2] and also in Fig. H] below. 

The evolution of a system towards a final equilibrium distribution as a function of po- 
tential well depth D is shown in Fig. [3l The parameters are similar to the Ar/W system of 
Fig. [1] above, the projectile is incident normally at 9i = 0° with Ei = 1 meV and Ts = 303 K 
but completely converged calculations (meaning less than 1% of the incident particles remain 
trapped) are shown for the three different well depths of 20 meV, 50 meV and 80 meV. This 
figure shows clearly that even if the incident energy is very small and the initially trapped 
fraction is large, the total scattered distribution does not approach equilibrium unless the 
adsorption well is sufficiently deep to cause long average trapping times. For a shallow well 
of 20 meV with an average desorption time of about 5.3 x 10~^° s as shown in Table [H The 
final distribution deviates strongly from an equilibrium distribution. It is only when the 
well depth is increased to about 80 meV, with the corresponding average trapping time of 
about 4.4 X 10~^ s, that near-equilibrium conditions are achieved. 

The response of the angular distribution of the final scattered particles for the same 
conditions as shown in Fig. [3] is given in Fig. HI The progression towards a cosine distribution 
is clearly evident with increasing well depth, but even for the quite large value of D = 200 
meV, for which the average trapping time is 6.4 x 10""^ s, the result deviates somewhat from 
cosine behavior. At Of = 0° the calculated value is 96% of the cosine maximum. 

An example showing the effect of mass on the convergence towards equilibrium is shown 
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in Fig. [5] corresponding to Ne scattering from tungsten at 1 meV of incident energy. For 
this system the mass ratio is roughly half that of Ar/W and approximate equilibrium be- 
havior is not achieved unless the well depth is approximately 150-200 meV in depth. The 
corresponding average trapping time for D = 200 meV is about 3.9 x 10^^ s as seen from 
Table [m 

The effects of mass on the approach to equilibrium is even more dramatically exhibited 
in Fig. [6] which shows the scattered energy distribution as a function of the projectile to 
surface mass ratio /i. The well depth is chosen to be 50 meV, the incident energy is 1 meV, 
6i = 45° and the temperature is 303 K. When the mass ratio is small, for instance fi = 0.005 
corresponds to hydrogen atoms scattering from tungsten, the total scattered intensity is far 
from an equilibrium distribution. However, for a mass ratio of about one-half the scattered 
spectra is very nearly in equilibrium. 

Larger temperatures tend to make initial trapping more difficult and lead to more rapid 
desorption, thus at high temperatures the trapping-desorption fraction will deviate more 
strongly from an equilibrium distribution. This is shown in Fig. [7] where calculations are 
shown for the same initial conditions as Fig. [1], i.e., low incident energy and a well depth 
of 80 meV, but three different surface temperatures of 303, 600 and 1200 K. For the lowest 
temperature the scattering is in very good agreement with the corresponding Knudsen equi- 
librium curve, but begins to deviate quite strongly as the temperature is increased. As the 
temperature is increased to the point where the trapping-desorption no longer is equilibrium, 
the average final energy becomes less than the 2kBTs value of the Knudsen distribution. 

The method of calculation presented here permits an examination of the energy distri- 
bution of the trapped particles at any average time after the initial collision. An example 
of this is shown in Fig. [8] for Ar/W with the same initial conditions as in Fig. [1] Both the 
negative energy trapped fraction and the positive energy chattering fraction are exhibited 
and it is seen that even after a very few iterations a smooth distribution with a maximum 
in the negative energy range appears. As the iteration number increases the most probable 
energy of all of the trapped particles shifts downward towards the bottom of the well as the 
trapped particles continue to lose energy to the surface on average. For very large numbers 
of iterations the trapped distribution begins to look exponential-like as a function of energy 
and reaches a steady state distribution that retains essentially the same functional form 
but gradually decreases in total integrated area as more and more particles are desorbed. 
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Interestingly, just after the initial collision the positive energy chattering fraction extends 
outwards to rather large energies with non-zero amplitude at positive total energies larger 
than the magnitude of the well depth. 

Some representative calculations of average trapping times r are shown in Tables HI and HTl 
based on the assumption of a square well of width 6 = 3 A. Table [I] shows calculations for 
the Ar/W system of Fig. [T]with 1 meV of incident energy and Ts = 303 K. The trapping 
times for several well depths are calculated in two slightly different ways, first by obtaining 
the average speed after each iteration as in Eqs. (fT9!) and fl20|) and then also from the root 
mean square speed. The two contributions to the total r from the negative energy trapped 
fraction and from the positive energy trapped fraction are exhibited separately. Also shown 
in Table [T] are the trapped fractions after the initial collision P°. Table [III shows similar 
information calculated for the Ne/W system of Fig. [5] with the same incident energy and 
temperature. 

For both systems the results are similar. For shallow well depths the trapping times 
are very short and the average time spent in negative and positive energy trapped states is 
comparable. As the well depth is increased, trapping times increase dramatically and the 
average time spent in the positive energy chattering state becomes negligible compared to 
the average time in the negative energy bound states. This increase in trapping times is 
nearly exponential as a function of well depth as is seen in Fig. [9] which graphs the numbers 
presented in Tables H] and [Tll Generally, the trapping times based on the rms speed after 
each iteration are somewhat smaller than those based on a calculation of the average speed. 
For Ar/W with a physically reasonable well depth of around 100 meV the average trapping 
time for 99% of the initially trapped particles to desorb is approximately 10~^ s. 

All of the above calculations have been for energies relatively small compared to the 
well depths in order to illustrate the conditions for which the trapping-desorption fraction 
approaches an equilibrium distribution. When the incident energy becomes comparable to 
or larger than the well depth the nature of the scattered intensity becomes quite different. 
As noticed in an important series of experiments first performed some years ago the intensity 
often exhibits a double-peaked structure, with a high energy peak due to direct scattering 
and a lower energy peak arising from the trapping-desorption fraction.^ Fig. [10] shows this 
for a system corresponding to Ar/W with 6i = 45°, D = 80 meV, Ts = 303 K and incident 
energies ranging from 100 to 500 meV as marked. The two contributions are shown in 
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separate panels, with the upper panel giving the direct scattering and the lower panel the 
trapping-desorption. With increasing Ei the total integrated direct scattering becomes larger 
and the peak becomes broader, the width roughly increasing proportionately to On 
the other hand the total trapping-desorption fraction becomes smaller and the shape of the 
distribution becomes less and less like that of an equilibrium Knudsen curve. For large Ei 
the trapping desorption intensity develops a long high-energy tail although its peak position 
always remains near to the most probable energy of the Knudsen distribution. The situation 
in which the differential reflection coefficient exhibits both distinct direct and a trapping- 
desorption peaks is discussed more in the next section in the context of comparisons of the 
present theoretical model with recent experimental data. 

IV. COMPARISONS WITH EXPERIMENT 

The calculations exhibited in the above section describe the range of initial conditions that 
lead to an equilibrium distribution in the trapping-desorption fraction but no comparisons 
with experimental data other than to a Knudsen distribution were made. However, in order 
to be credible, the theoretical approach should be capable of explaining real experiments. 
Demonstrating that ability is the objective of this section. 

There is a long history of gas-surface scattering experiments using hyperthermal 
atoms^^^^^^^-^^^^^^^^ and moleculea^l^^^^^^^^^^^^^ as projectiles. 
If the projectile gas has mass larger than that of hydrogen or helium such high energies im- 
ply that the scattering will be classical, which means that many phonons will be transferred 
in the collision. This is the type of experiment that should be amenable to the theoretical 
treatment described here. In many cases the energy-resolved scattered spectra exhibit a 
double peaked structure, with a somewhat narrow high-energy peak centered at smaller en- 
ergy than the incident beam energy (if the incident energy is large compared to the surface 
temperature) and a broader low-energy peak at thermal energies. The usual interpretation 
has been that the high-energy peak is direct scattering from a single collision (or at most, 
a very small number of collisions) and the low-energy peak arises from trapping in the ph- 
ysisorption well of the interaction potential with subsequent desorption at a sufficiently later 
time so that those particles come into near equilibrium at the surface temperature.- 

A recent and important paper reporting extensive measurements that show clearly a 



15 



set of conditions for which direct scattering and trapping-desorption can be observed is 
that of Gibson, Isa and Sibener for scattering of Ar from an ordered 1-decanethiol self- 
assembled overlayer on a Au(lll) substrate.- The experiments were carried out with well- 
defined monoenergetic beams of Ar incident at energies ranging from roughly 60 to 600 
meV, and with both incident and detector angles independently variable and ranging from 
near-normal to near-grazing with respect to the surface. All measurements were made in 
the scattering plane (the plane containing the surface normal and the incident beam) which 
was aligned along the (110) direction of the Au(lll) surface. At low incident energies 
and if 9i or Of was near-normal they did not observe a clear double peaked intensity in 
the scattered spectra. However, at higher energies and for large incident or final angles 
the characteristic double-peaked structure was very apparent. They analyzed their data 
quantitatively with an ad hoc model consisting of the sum of a shifted Maxwell-Boltzmann 
distribution to fit the direct scattering and an equilibrium distribution to fit the trapping- 
desorption fraction. They also made some more qualitative analysis of their data using 
classical trajectory calculations developed by Hase and coworkers.— In the process of their 
analysis they determined, by assuming that the direct scattering was due to a single collision 
and using well-known Baule relations for binary collisions, that the effective mass of the 
surface imphed a mass ratio n = 0.62, or Mc = 64.4 amu as opposed to the total mass of 
the 1-decanethiol which is 174.3 amu. The potential energy landscape function developed 
for the classical trajectory calculations had a well depth ranging from 33 meV at the on-top 
sites above the terminal CH3 groups to 67 meV in the center of the rhombus formed by a 
group of four of the methyl groups.- 

Fig. [11] shows an example of calculations compared to the Ar scattering data taken from 
the upper panel of Fig. 2 in Ref. [^] and is for their lowest incident energy Ei = 65.3 meV. 
The data was reported as intensity versus time-of-fiight (TOF) and the calculations have 
been transformed accordingly. The other incident parameters are 6i = 45°, Of = 50° and 
Ts = 135 K. The calculations were carried out for a well depth D = 35 meV and an effective 
surface mass Mc = 71 amu {fi = 0.56). 

This effective surface mass ratio is slightly smaller than the value fi = 0.62 estimated in 
Ref. [^]. For this case of low incident energy the calculations are not particularly sensitive 
to the value of fi because there is very little evidence for a significant direct scattering 
component. However the value of /i = 0.56 is chosen as a consequence of comparisons with 
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the higher energy data shown below, where the much more pronounced direct scattering 
component is extremely sensitive to fi. In this example the data do not exhibit a double- 
peaked structure and the calculated most probable energy (peak position) has a TOF time 
corresponding to Ef = 12.57 meV, very close to the equilibrium value of 11.6 meV at this 
temperature, again indicating that the scattering is mostly trapping-desorption. 

Three calculated curves are shown in Fig. [11], the solid curve is the total differential reflec- 
tion coefficient converted to TOF, the dashed curve is the trapping-desorption contribution 
only, and for comparison a Knudsen equilibrium distribution is included as a dash-dot curve. 
The Knudsen distribution and the calculation are both normalized to unit intensity as in 
Eq. ([T]). Although it is apparent that the Knudsen curve, if renormalized to fit as closely as 
possible to the data, would match essentially as well as the total calculated intensity, the 
fact that it is smaller and nearly the same as the calculated trapping-desorption fraction 
indicates that there is significant direct scattering but its most probable energy and width 
is nearly the same as the equilibrium distribution. Because of the strong overlap of the 
direct and trapping-desorption fractions, it is not surprising that the data of Fig. [H] can be 
matched roughly as well by a total scattering intensity using a range of well depths from 
20 to somewhat over 35 meV. We have chosen D = 35 meV because of the much stronger 
constraints placed on this parameter by the higher energy data considered below. 

Three examples of data measured at the intermediate energy of 365 meV are shown in 
Fig. [12] at the same temperature of 135 K and for three different combinations of incident and 
final angles. The middle panel for Of = 50° and 6i = 30°, relatively close to normal incidence, 
does not exhibit a double peaked structure. The other two panels, for 6i = 45° and Of = 50° 
(upper panel) and 6i = 30° and Of = 80° (lower panel) present a clear distinction between 
the rather sharp peak at short TOF and a broader shoulder at larger times. These data 
were taken from the middle panel of Fig. 2 and the lower two panels of Fig. 4, respectively, 
of Ref. p]. The solid curves in Fig. [12] are calculations carried out with fi = 0.56 and = 35 
meV. The calculations explain the data quite well, and they show clearly the separation 
between the direct and trapping-desorption fractions. The value D = 35 meV agrees well 
with the that of the potential energy function for this system developed in Ref. p] 

Also shown in Fig.[T2]are the trapping-desorption fraction and the Knudsen curves. Inter- 
estingly, the trapping-desorption fraction itself has a multiple-peaked structure with a small 
sub-peak appearing at almost the same final energy as the direct scattering contribution. 
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This small high-energy sub-peak comes from the first few collisions as the initially adsorbed 
particles travel in the potential well. These first few collisions have a high probability of 
ejecting particles back into the continuum with relatively little loss of energy compared to 
the direct scattering fraction. However, it is clear that the largest part of the trapping- 
desorption fraction resembles closely the shape of the Knudsen curve. 

It also becomes apparent from Fig. [12] that there is a straightforward manner in which the 
comparison of calculations with data allows for the determination of the two parameters. The 
effective mass determines the most probable final energy of the direct contribution and then 
the well depth determines the relative intensity of the trapping-desorption fraction which 
becomes bigger with increasing D. The peak position of the direct scattering contribution is 
extremely sensitive to the mass ratio and this is why we chose the value /i = 0.56 as opposed 
to the value ^ = 0.62 of Ref. ^] which was based on the Baule equations describing hard 
sphere scattering. 

Three examples of data for scattering at the high energy Ei = 582 meV, all of which 
exhibit the double-peaked structure, are shown in Fig. [131 The data were taken from the 
lower panel of Fig. 2 and the middle and lower panel of Fig. 3, respectively, of Ref. [^]. The 
upper panel of Fig. [13] is for 6i = 45° and Of = 50°, the middle panel is for 6i = 45° and 
6*/ = 40°, and the lower panel is for 9i = 60° and 9f = 40°. 

In the upper two panels with 6i = 45° two curves showing the total scattering intensity 
are shown, for D = 35 and 45 meV while in the lower panel with 6i = 60° only tht = 35 
meV calculation is shown. All calculations were done with /x = 0.56, the value that leads to 
agreement with the data for the direct scattering peak. The dotted curve is the Knudsen 
distribution and the dash-dotted curves show the trapping-desorption fraction. 

It is interesting that at this larger incident energy the two cases with the more normal 
incident angle of 45° require a well depth of 45 meV in order to obtain agreement between 
calculations and data, while for the much more grazing incidence of 80° the best well depth 
is 35 meV, the same as used for all the lower energy calculations. This appears to indicate 
that for larger normal incident energy the incoming atoms are probing deeper parts of the 
potential energy landscape, which as mentioned above has been estimated to have a well 
with depths that vary between 33 and 67 meV.- Again, as in Fig. [12] the trapping-desorption 
fraction exhibits structure at larger final energies near the energy of the direct scattering. 

Comparison of the present calculations to this Ar scattering data leads to a few general 
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comments that can be applied to the observed energy-resolved spectra for cases in which a 
double structure appears due to the distinct phenomena of direct scattering and trapping- 
desorption: 

(1) A characteristic double-feature structure, with a well-defined direct scattering peak and a 
secondary peak or shoulder arising from trapping-desorption appears only at relatively high 
incident beam energy and when one or the other of 9i or 6'/ is large, as was already clear from 
Ref. [^]. It is also necessary that the physisorption well depth is sufficiently large to cause 
significant trapping in the bound states during the initial collision, and the temperature 
must be smaller than D. However, for near normal incidence conditions and with a large 
well, the initial trapping becomes so large that the direct scattering contribution becomes 
small, and this explains the need for the incident angle to be relatively large. This situation 
becomes evident in Fig. [T3] which shows calculations for Ei = 582 meV with 6i = Of = 10° 
and -D = 35 meV. Even though the incident energy is large compared to the well depth, the 
normal incidence conditions gives rise to such large trapping that the total scattering is not 
very different from the trapping-desorption fraction. 

(2) When the direct scattering contribution is significant, the trapping-desorption inten- 
sity deviates substantially from that of an equilibrium Knudsen distribution. In fact, the 
trapping-desorption signal can exhibit structure and small peak-like features at high energies 
close to those of the direct contribution as shown in Figs. [12] and [131 However, in many cases 
the largest portion of the trapping-desorption intensity resembles the shape of a Knudsen 
distribution and this is especially true for the low energy tail. 

(3) It is interesting to note that the direct and trapping-desorption fractions should have 
very characteristic and quite different signature behaviors in their temperature dependence. 
The direct scattering, which in the present calculations arises from a single collision with 
the surface, is essentially given by Eq. ([T8|) . Under conditions for which the direct scattering 
of Eq. f[T8|) appears nearly Gaussian-like in the energy transfer, which is the situation in 
several of the cases shown here, the mean square energy deviation which is proportional to 
the full width at half-maximum (FWHM) isl^diiii. 

= = ^^^^^^^^ ' ^^^^ 
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where 



1^1 COS 9 

with 



where 6^ is the total scattering angle (i.e., the angle between pj and p/) and Eqs. fl23|) - fl25|) are 
obtained under the assumption of binary collision conditions for which Ef = f{fi, 9)Ei. Thus, 
Eq. fl23|) shows that the FWHM of the direct scattering peak will be approximately propor- 
tional to v^Ts, which is the characteristic of the multiphonon scattering regime. However, 
the trapping-desorption fraction will have a FWHM temperature dependence more closely 
approximating the linear in T5 behavior of the Knudsen distribution of Eq. fl22l) . Similarly, 
the most probable intensity (maximum peak intensity) of the direct scattering will vary as 
Ijy/Ts according to Eq. ffTSl) while the trapping-desorption peak intensity should behave 
more like that of the I/T5 behavior of the Knudsen distribution of Eq. (!22l) . 
(4) Finally, the comparison of the present calculations with the data provide a simple way 
to extract the two relevant parameters. The position in final energy of the narrow direct 
scattering peak determines the value of the effective surface mass, and in fact, the most 
probable energy is quite sensitive to this parameter. The intensity of the broad trapping- 
desorption peak increases with the well depth, and fixing the relative intensities of the two 
contributions determines D. This indicates that for the large incident energies considered 
here the primary influence of the well depth is to establish the initial trapping fraction. Once 
trapped, the details of the shape of the potential well are not important as is evidenced by 
the fact that the long-time tail trapping-desorption fraction (the low energy tail) eventually 
desorbs at thermal energies. 



V. CONCLUSIONS 



In this paper we have developed a theoretical formulation of the scattering of atomic 
projectiles with surfaces that includes not only the direct scattering arising from a single, 
or a small number of collisions with the surface but also allows for trapping and subsequent 
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collisions of trapped particles inside the physisorption well. The trapped particles can be 
followed until they eventually all desorb and leave the surface region. The multiple coUisions 
of the initially trapped fraction with the surface are treated with an iteration algorithm that 
tracks trapped particles with both negative and positive total energies and determines at 
each subsequent collision the fraction scattered back into the positive energy continuum 
which then leaves the surface region. 

Using this theoretical formalism we have calculated numerous examples, firstly in order 
to establish the conditions under which the trapped and subsequently desorbed particles ap- 
proach an equilibrium distribution, i.e., to establish the conditions under which the Maxwell 
assumption is valid, and secondly wc have used the theory to produce quantitative agreement 
with recent measurements thus providing explanations of the basic underlying processes that 
give rise to the experimental scattered spectra. 

Under many conditions, the observed spectra in gas-surface scattering experiments consist 
of two distinct contributions. The first contribution is the direct scattering part which 
is usually a relatively sharp peak with a most probable energy somewhat lower than the 
incident energy provided that the surface temperature is not large compared to the incident 
energy. The second of these contributions is the trapping-desorption, attributed to particles 
that are initially trapped and then spend a large time moving in the physisorption well 
where they slowly begin to exchange energy with the surface and then eventually desorb in 
a distribution at thermal energies corresponding to the surface temperature. A large part of 
the work considered here is devoted to determining when the trapping-desorption fraction 
approaches an equilibrium Knudsen flux. 

Our calculations show that under conditions in which a clear direct scattering and 
trapping-desorption double-peaked structure is evident in the energy-resolved spectra the 
trapping-desorption fraction, although mainly emitted at thermal energies, can differ con- 
siderably from an equilibrium distribution. It can even exhibit structure consisting of small 
peaks at high energy near the most probable energy of the direct scattered intensity. These 
higher energy peaks arise from the first few collisions with the surface inside the well in- 
dicating that these initial collisions have a large probability of ejecting particles into the 
continuum with little additional energy loss as compared to the direct scattering. However, 
even under conditions for which the trapping-desorption fraction is highly non-equilibrium 
its low energy tail still is well described by a Knudsen distribution. Thus, our calculations 
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show that the Maxwell assumption is rarely achieved in real experimental conditions, al- 
though it is very useful as an approximate guide as is evident from the fact that it is still 
often used as a method to analyze measured data. 

We carried out a number of calculations in order to characterize the conditions under 
which the trapping-desorption fraction does approach an equilibrium distribution. Basically, 
equilibrium behavior is achieved only for cases where the direct scattering is negligible 
and average trapping times are long, which implies that nearly all of the incident beam is 
adsorbed after the first collision. This implies an incident energy relatively small compared 
to the interaction potential physisorption well depth and temperatures corresponding to 
energies (measured in units of ksTs) that are also small compared to the well depth. The 
approach to equilibrium occurs more rapidly with larger gas-to-surface-atom mass ratios 
when this ratio is less than unity. 

The approach to equilibrium of the trapping-desorption fraction was studied as a func- 
tion of of all the initial experimental parameters that can be manipulated, including the 
projectile and surface mass, the well depth of the potential, the incident energy and angles, 
the final scattering angles and the surface temperature. For example, as the well depth is 
increased with all other parameters held constant, wc find that the energy-resolved scattered 
spectrum rather quickly approaches that of a Knudsen distribution when the well depth be- 
comes significantly larger than the incident energy provided the temperature is also small 
compared to the well depth. The angular behavior becomes independent of azimuthal angle 
under the same conditions that the energy distribution becomes Knudsen-like. However, 
the polar angle cos Of behavior is only approximately obeyed for conditions under which 
the energy dependence first approaches equilibrium, even though the energy dependence is 
nearly Knudsen-like at all polar angles. Only for well depths very large compared to the 
incident energy and temperature does the trapping-desorption fraction achieve the classic 
Knudsen cos9f shape. 

In this formalism, because of its iterative approach, the number of collisions of the trapped 
particles with the surface can be followed. This means that trapping times can be calculated 
as well as other information such as the relative fractions of particles trapped with negative 
total energies and those trapped in the chattering states having positive total energies. 
Under conditions in which the trapping-desorption fraction did achieve near equilibrium the 
trapping times were estimated to be as large as 10~^ — 10"^ s. The positive total energy 
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chattering fraction can be large under conditions where equihbrium is not achieved, but 
for conditions that produce equihbrium in the trapping-desorption the fraction of trapped 
particles residing in the chattering states becomes negligible. 

An important aspect of this work is that the theoretical model used for the calculations 
provides quantitative explanations of experimental measurements. We made comparisons 
with important and recent Ar scattering data obtained in beams-surface scattering exper- 
iments with a surface consisting of a self-assembled adsorbed layer of 1-decanethiol on a 
well-ordered Au(lll) substrate.^ In agreement with the experimental observations, we found 
that clearly distinguishable direct and trapping-desorption contributions arise when the in- 
cident beam energy is large compared to the potential well depth and when one or the other 
of the incident or detector angles was large relative to the surface normal. In addition the 
calculations indicate that, in order to resolve distinct direct and trapping-desorption fea- 
tures, the temperature must be smah compared to the well depth and the well depth must 
be large enough to cause trapping of a significant fraction of the incident beam at the initial 
collision. 

Under conditions for which distinct direct and trapping-desorption features were evident, 
the data can be used to determine two important characteristics of the interaction potential, 
the effective surface mass of the adsorbate and the well depth. The effective mass is deter- 
mined by matching the calculated direct scattering peak to that of the experiment, and the 
well depth then is determined by matching the relative intensity of the trapping-desorption 
contribution. 

An interesting prediction coming out of this work is that the direct and trapping- 
desorption contributions have clearly different signature behaviors as a function of surface 
temperature. The FWHM of the direct peak should increase approximately with the square 
root of the temperature whereas the trapping-desorption has a full width that increases 
approximately linearly with Ts, similarly to the Knudsen distribution. The most probable 
intensity of the direct peak, according to the scattering model used here, decreases inversely 
with the square root of temperature, while the trapping-desorption decreases roughly lin- 
early with the inverse of the temperature, again similarly to the Knudsen distribution. For 
both peaks, the increase in FWHM and decrease of most probable intensity is the behavior 
expected in order to preserve unitarity. It should be noted that the direct scattering, as 
shown in the approximations to Eq. ( JTSl) appearing in Eqs. (1231) -( l25l) . exhibits the same 
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behavior in the incident energy Ei as it does in the temperature T5. However, the most 
appropriate parameter in which to carry out an experimental search for these behaviors 
would be the temperature, because the interaction potential surface is likely to change with 
variation of Ei but is less likely to change with T5. 

This work demonstrates that calculations of direct scattering, trapping and desorption in 
atom-surface scattering can provide real quantitative explanations of experiments as well as 
indicating the conditions for the validity of Maxwell's assumption on the equilibrium nature 
of the trapping-desorption fraction. It shows that the interaction potential model must 
contain two essential ingredients: a physisorption well depth and it must allow for transfer 
of mechanical energy between the projectile and the surface atoms. However, it also shows 
that the most important aspect is to have a theory in which the statistical mechanics is 
treated in a reasonably correct manner. 
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TABLE L The desorption times r and initial sticking fractions P° for Ar/W. Ei is 1 meV, 9i is 
45°, Ts is 303K and well depths ranging from 20 meV to 200 meV are shown. The upper set of 
values was obtained from the normal rms speed of trapped particles and the lower set are values 
determined from the average normal speed. 
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TABLE II: The desorption time for Ne/W with other parameters the same as in Table HI 
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FIG. 1: Argon scattering from a tungsten surface: the evolution of the final energy distribution 
for particles scattered into the continuum states after a specified number of iterations. The surface 
temperature is 303 K, the incident energy Ei = 1 meV, the well depth D is 80 meV and the incident 
angle 0i = 45°. Five curves of the final distributions for the iteration numbers iV = 1, 5, 50, 500 
and 2124 are shown. For comparison a Knudsen equilibrium flux is given as open circles. 
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FIG. 2: The evolution of the final distribution in polar angle 9f for the Ar/W system with the 
same parameters as in Fig. [TJ A Knudsen cosine distribution is shown as open circles. 
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FIG. 3: Final energy distributions Ar/W system as a function of potential well depths D = 20 
meV, 50 meV and 80 meV as shown. The temperature of the surface is 303 K, the incident energy 
Ei = 1 meV and 6i = 0°. A Knudsen distribution is shown as open circles. 
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FIG. 4: The polar angular distribution for the same system shown in Fig. [3l The evolution towards 
the Knudsen distribution, displayed as open circles, is shown for a series of increasing well depths. 
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FIG. 5: Energy distribution of the scattered particle for Ne/W for several well depths as shown. 
The temperature of the surface is 303 K, the incident energy Ei = 1 meV and 6i = 0°. The 
Knudsen distribution is shown as open circles. 
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FIG. 6: The energy distribution as a function of the mass ratio /x for a system with Oi = 45°, 
Ei = 1 meV, Ts = 303 K and i:> = 50 meV. 
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FIG. 7: Energy distribution of the scattered particles for Ar/W with incident energy Ei = 1 meV, 
6i = 45° and a well depth D = 80meV. Three different surface temperatures Ts = 303 K, 600 
K and 1200 K are shown as curves and the corresponding equilibrium distributions are shown as 
points. 
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FIG. 8: The evolution of the energy distribution inside the potential well as a function of iteration 
number for the same system as in Fig. [TJ 
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FIG. 9: The desorption times for Ar/W and Ne/W as functions of the well depth for the incident 
conditions shown in Tables HI and ITTl 
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FIG. 10: The direct and trapping-desorption scattering energy distributions for incident energies 
large compared to the well depth. The system is Ar/W with Oi = 45°, D = 80 meV, T5 = 303 
K and incident energies as marked, a) the upper panel shows the direct scattering contribution, 
and b) the lower panel shows the trapping-desorption intensity. A Knudsen distribution is also 
included in the lower panel. 
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FIG. 11: Intensity versus TOF for Ar scattering from a 1-decanethiol layer on Au(lll) with 
Ei = 65.3 meV, 6i = 45° and 6f = 50°. The calculation shown as a solid curve is the total 
differential reflection coefficient converted to TOF calculated with fi = 0.56 and D = 35 meV, 
the dash-dotted curve is the trapping-desorption fraction, and the dotted curve is a Knudsen 
equilibrium distribution. 
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FIG. 12: Intensity versus TOF for Ar scattering from a 1-decanethiol layer on Au(lll) with 
Ei = 356 meV: a) 9i = 45° and 9f = 50°, b) Oi = 30° and 0f = 50° and c) Oi = 30° and 0/ = 80°. 
The solid curves are calculations with fi = 0.56 and D = 35 meV, the dash-dotted curves are the 
trapping-desorption fractions and the dotted curves are the Knudsen distribution. 
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FIG. 13: Intensity versus TOF for Ar scattering from a 1-decanethiol layer on Au(lll) with 
Ei = 582 meV: a) Oi = 45° and Of = 50°, b) 9i = 45° and 9f = 40° and c) 9i = 60° and (9/ = 40°. 
The curves are labeled the same as in Fig. [12] except calculations are shown for both D = 35 and 
45 meV. 
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FIG. 14: Calculations similar to those of Fig. [13] with Ei = 582 meV and D = 35 meV but with 
6i = 10° and Of = 10°, showing that for normal incidence a distinct direct scattering peak is not 
expected. 
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